Interaction-disorder competition in a spin system evaluated through the Loschmidt 
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The interplay between interactions and disorder in closed quantum many-body systems is rel- 
evant for thermalization phenomenon. In this article, we address this competition in an infinite 
temperature spin system, by means of the Loschmidt Echo (LE) , which is based on a time reversal 
procedure. This quantity has been formerly employed to connect quantum and classical chaos, and 
in the present many-body scenario we use it as a dynamical witness. We assess the LE time scales as 
function of disorder and interaction strenghs. The strategy enables a qualitative phase diagram that 
shows the regions of ergodic and non-ergodic behavior of the polarization under the echo dynamics. 
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The ergodic hypothesis of statistical mechanics implies 
the equivalence between time and ensemble averages. It 
is expected that a conservative many-body system satis- 
fying such hypothesis would explore uniformly the whole 
energy shell. It is now a long time since Fermi, Pasta 
and Ulam (FPU) questioned how the irregular dy- 
namics induced by non-linearities in a Hamiltonian may 
lead to energy equipartition as a manifestation of ergod- 
icity. Such dynamics, referred as thermalization, did not 
show up in their pioneering numerical simulations. Even 
though their striking results are now explained by the 
theory of chaos 0] , the solution of the quantum analogues 
are still in early stages 0, Q • 

Thermalization in isolated, strongly interacting, quan- 
tum systems is defined relative to a certain set of ob- 
servables [^. In particular, remarkable experiments have 
been performed to monitor momentum distributions of 
cold atoms loaded in one-dimensional optical lattices , 
where the integrability of the underlying dynamics is 
weakly broken. Accordingly, a fundamental question is 
whether a non-ergodic to ergodic transition threshold ex- 
ists as one may go parametrically from an integrable to a 
non-integrable quantum system. In the FPU problem one 
can surely answer affirmatively, since the onset of dynam- 
ical chaos can play the role for such transition [2|. For 
interacting quantum systems in the presence of disorder, 
it is expected that a phase transition exists between de- 
localized states and a phase characterized by many-body 
localization (MBL) Q. This would constitute the sought 
threshold between ergodic and non-ergodic behavior 

The MBL results from a quantum dynamical phase 
transition between diffusive (ergodic) and localized (non- 
ergodic) dynamics 0] ■ It occurs at nonzero (and even- 
tually infinite) temperature, and is manifested in dynam- 
ical properties. The crucial idea is a competition be- 
tween interactions and Anderson disorder [lO|. In the 
many-body ergodic phase, the expectation values of ob- 
servables computed on a finite subsystem, using a single 



energy eigenstate of the whole interacting system, would 
coincide with those evaluated in the corresponding mi- 
crocanonical thermal ensemble UJ, |l2| . In this condition 
one may say that the system acts as its own heat bath. 
Quite on the contrary, in the non-ergodic many-body 
phase, dynamics resembles a "glassy" behaviour that pre- 
cludes self-thermalization. In view of the difficulties of 
addressing a full many-body dynamics of specific corre- 



lation functions 13|, , much of the progress in assess- 



ing the transition between the mentioned regimes relied 
mainly on the evaluation of spectral properties • 

The dynamics of specific observables has probed useful 
to monitor the onset of many-body chaos [l6| . In par- 
ticular, the Survival Probability (SP) of the eigenstates 
of an unperturbed Hamiltonian under the evolution 
of the full interacting Hamiltonian Hq + S decays with a 
characteristic time scale r. In such case, a crossover from 
an exponential controlled by l/r cx to a gaussian 

characterized by l/r cx ||S|| is interpreted as evidence 
of the onset of a chaotic structure in the eigenstates of 
Ha + Yi [rzl - flit . The first regime is described by the 
Fermi golden rule (FGR), and hence the breakdown of 
the FGR leads to dynamical chaos. Much in the spirit 
of these dynamical approach, in this Letter we propose 
and analyze the evolution of an experimentally accessi- 
ble local observable, as a way to assess the interaction- 
disorder competition in a many- spin system. We resort 
to the Loschmidt Echo (LE) [20| , a measure of the revival 
that occurs when a time-reversal procedure is applied to 
Hq. The LE has been experimentally evaluated from the 
local polarization of spin systems to quantify the role of 
perturbations (i.e. E) and the system's own complex- 



ity on dynamical reversibility [21[. Besides, in classically 



chaotic systems, it is well-known that the LE of a semi- 
classical excitation undergoes a transition into a regime 
where its decay rate is given by the classical Lyapunov 
exponent [13]. Here, we use the LE dynamics as a time- 
dependent autocorrelation function that quantifies the 



2 



time-scales and ergodicity. Specifically, the time reversal 
procedure "filters out" the irrelevant dynamics produced 
by an integrable Hq which would hide the sought informa- 
tion [2^. Thus, the LE becomes a privileged dynamical 
witness of the competition between the interactions and 
disorder and a potential experimental candidate. 

The spin model- We consider a spin chain, whose dy- 
namics is given by the total Hamiltonian H — Ho + T,: 



N 

Ho=Y.j[S^Sf+,+SfSf^,] 
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where A is the magnitude of the homogeneous Ising 
interaction and hi are randomly distributed fields in an 
interval [—h,h]. Periodic boundary conditions (ring) are 
imposed and, unless explicitly stated, N — 12. No- 
tice that Hq can be mapped into two independent non- 
interacting fermion systems by the Wigner- Jordan trans- 
formation 24] , while S includes both the two-body Ising 
interaction and the local fields (disorder). Since Hq en- 
closes single-particle physics, we consider it as the irrele- 
vant part of the dynamics, a term that the LE manages 
to get rid of, allowing to focus on interactions and disor- 
der. This fact justifies the idea of using the LE as a filter 
for the relevant physical processes. 

The spin autocorrelation function evaluated as a (local) 
Loschmidt Echo.- We consider a high (infinite) temper- 
ature state formally denoted by l^'eg), which represents 
an ensemble average over all basis states with the same 
statistical weights. We study the spin autocorrelation 
function at the particular site 1, 



Mi,i(2tfl) 



{^,g\Sf{2tn)Sm\^eg) 



(3) 



where the spin operator is written in the Heisenberg pic- 
ture as: 



S',{2tR) = ul{tR)Ul{tR)SfU_{ti,)U+{tR). (4) 

The evolution operators are U^{tR) = cxp[— -^(//q + 
f:)tR] and UjtR) = exp[-i(-i?o + ±)tR]. There- 
fore, it is explicit that the echo procedure performed 
over Hq yields a global evolution operator U{2tR) — 
U^{tR)U^{tR). Using cyclic invariance of the trace, one 
can replace the average over all basis states by an av- 
erage over all states that have spin 1 up-polarized [23| . 
Additionally, we replace the ensemble average by an av- 
erage over a few entangled states [2^, which provides 
a quadratic speedup of computational efforts. It yields 



FIG. 1. Color online. Loschmidt Echo for different parameter 
regimes. The characteristic time scale is defined by the decay 
up to 2/3, shown by the horizontal dotted line, (a) ft = 0, 
< A < 1.5J. Smooth decay, until saturation is reached 
(horizontal dashed line), (b) h ^ 0, 1.5J < A < 4.8J. LE 
long tails destroys saturation. The SP given by Eq. [7| de- 
termines the limit time scale (shaded region), (c) A = 0, 
< h < 1.0 J. The localization length remains bigger than 
the system size, (d) A = 0, l.OJ < h < 5.5J. The local- 
ization length turns to be smaller than the system's size, and 
the polarization keeps around site 1. 



equivalent results provided that only local observables 
(e.g. Eq. [3]) are evaluated. Thus, we consider: 



\^ne,) = Itl) 



JV-1 



(5) 



where t/?^ is a random phase and {|/3r)} are state vectors 
in the computational Ising basis of the — 1 remaining 
spins. The Mi 1 is now written in the Schrodinger picture 
as: 



Mi,i(2t) 



2 {^ne,\ Ul{t)Ul{t)SlU_{t)U4t) |*„e,/ 



(6) 

The explicit time dependence of Mi_i(2t) is evaluated 
by means of a 4th order Trotter-Suzuki decomposition 
without any Hilbert space truncation, implemented on 
General Purpose Graphical Processing Units f26| . 

Results - We start the evaluation of Eq. [S] with h = 
(no disorder), increasing the interaction strength A from 
zero, see Fig. [T] (a) and (b). The LE short-time scaling 
is 1 — Mi_i(t) oc t*, which can be analytically verified 
by expanding the evolution operators up to 4*'' order. 
Beyond the short-time window, it has a smooth decay 
produced by the non-reversed terms E [23| . As we do not 
have an explicit functional form for this regime, we define 
an effective characteristic time r for the plotted curves 
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as the decay up to 2/3. The rates 1/t are plotted in Fig. 
[2] as function of A, for different disorder magnitudes h. 

Beyond the decay regime, as shown in Fig. [l]-(a), 
the LE saturates at 1/N, which means that the initial 
polarization (local excitation defined by Eq. [5]) is uni- 
formly spread over the whole spin set. This is indeed 
the standard picture of decoherence process leading to 
an irreversible spread. But, if A is further increased, the 
LE-decay slows down showing long tails, see Fig. [l]-(b). 
These destroy the saturation at least in the time- window 
analyzed in the present work. Such regime may be as- 
sociated with a glassy polarization dynamics, i.e. the 
prevalence of the freezing effect of the Ising terms over 
the spreading induced by the XY ones. 

If we consider A = and let the disorder h increase, 
the picture is indeed the standard Anderson Localization 
problem. In such case, the localization length must be 
compared to the finite size of the system. Hence, for very 
weak disorder, the LE degrades smoothly, as the local- 
ization length is longer than the system's size. When 
the disorder is strong enough, the localization length is 
smaller than the system's size and thus the initial lo- 
cal excitation remains around the site 1. In fact, the 
crossover between these two physical situations can be 
quantified equating the system's size with the localiza- 
tion length A ~ 24J^//i^, given in a FGR estimation 



27|. This yields h ~ \/2J , in fairly good agreement with 
the behavior shown in Figs. [T] (c) and (d). 

Notice that in any case, neither Ising interactions nor 
Anderson disorder, can produce a LE-decay faster than a 
well defined time scale (see Figs. [11(b) and (d)). This is 
specifically determined by the SP of the local excitation 
under the evolution given by Hq, 

Pl^{2t) = 2 (*„e,| exp {im/h) SI exp [-im/h) |*„eg) 

(7) 

where we emphasize that there is no dependence on E. 
Such iJg-controUed decay resembles the perturbation in- 
dependent decay experimentally observed in spin systems 
[2lj and the Lyapunov regime of classically chaotic sys- 
tems (2^ . 

In Fig. [5]we show the role of the time-scale determined 
by Eq. [71 acting as the limit for LE decay rates. These 
results are quite general, similar behavior is obtained for 
next nearest neighbors interactions (both in _ffo and E), 
provided that E has only Ising terms or Anderson disor- 
der. 

In order to analyze the ergodicity of the polarization 
dynamics observed in our finite system, we evaluate the 
mean LE, MiX- 



Mi,i(r) = - / Mi,i{t)dt. 

^ JQ 

The standard analysis of the Anderson localization 
problem should imply, in present case, to compute 
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FIG. 2. Color online. LE decay rates as a function of the 
interaction strength A, for different disorder magnitudes h. 
Horizontal dashed line represents the (-ffo) SP time scale, i.e. 
Eq. [71 Data for A = was obtained from single-particle 
physics and 500 disorder realizations, while data for A > 
was obtained from full many-body simulations and 10 disor- 
der realizations. The smoothness of matching evidences the 
robustness of the simulations performed. 



liniT^oo Mi^i(T). Instead, we analyze Mi.i(T) at T = 
12h/J which, as shown in[Tl is long enough to allow for 
a uniform spreading of the polarization, provided that 
A is strong. Also, a rigorous upper bound for the in- 
tegration time T must be the system's Heisenberg time 
Tff, at which finite-size recurrences show up. In Fig. 
[3l we show a level plot of Mi,i as a function of the in- 
teraction A and disorder strength h. This results in a 
qualitative phase-diagram which evidences the competi- 
tion between such physical magnitudes. When both A 
and h are weak, Mi^i remains near 1, since the system is 
almost reversible. Thus, the reddish parametric region at 
the bottom left corner may be associated with decoher- 
ence, i.e. the system is weakly perturbed by uncontrolled 
degrees of freedom. If either A or /i are further increased, 
the system enters in a diffusive regime where the initial 
local polarization rapidly spreads irreversibly all across 
the spin system. Consistently, this bluish region is as- 
sociated with an ergodic behaviour for the polarization, 
since it is equally distributed along the spin system. For 
low disorder (h < 0.5 J), increasing A leads to an Ising 
predominance, which freezes the polarization dynamics. 
We interpret such behavior as a glassy dynamics, with 
long relaxation times. This localization keeps Mi^i high. 
Analogously, increasing h for a fixed value of A evidences 
a smooth crossover to a localized phase (green color), 
where the polarization does not diffuse considerably. 

Notice that localization by disorder is weakened when 
l.OJ < A < 2.0J, since the ergodic region seems to 
unfold for larger h. In fact, such values of interaction 
strength correspond to a faster arrival to the 1/iV satu- 
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ration, as shown in Figs. [T] (a) and (b) (when h = 0). 
Additionally, data around the A axis show that the disor- 
der tends to abruptly destroy the quenching produced by 
Ising interactions. This seems suggestive of a parameter 
region where the interaction-disorder competition leads 
to a sharp transition between glassy and ergodic phases. 
However, a reliable finite size scaling of this regime would 
require excesively long times to capture how a vitreous 
dynamics is affected by disorder. 

Previous numerical results of the SP 28[ and an anal- 
ysis of the many-body eigenstates Q of the same spin 
model have identified critical values for the MBL transi- 
tion. Quite remarkably, they lie precisely at the crossover 
between the ergodic and the localized phases of the LE; 
see data-points (a) and (b) in Fig. [3l In our simmula- 
tions, increasing N (e.g. 10, 12 and 14) enables a larger 
integration time T, since T < Th oc N. In fact, when 
A ^ l.OJ, it can be verified that both sides of the tran- 
sition are well behaved since Mi,i ^ 1/N in the ergodic 
regime, while Mi^i ~ 1/A for h strong enough (regardless 
of N), and dMi^i/dh increases with N . However, this fi- 
nite size scaling of Mi^i(r) within our accessible range 
is not enough to yield precise critical values for the MBL 
transition. 

In summary, we were able to draw a qualitative phase- 
diagram for the polarization under the LE dynamics, 
identifying ergodic, localized and glassy regimes. It dis- 
plays a non-trivial geography with a deep penetration 
of the ergodic phase into the glassy domain separating it 
from the localized region. Besides, while in finite 1-d sys- 
tems Anderson localization is indeed a smooth crossover, 
it seems to develop into a ergodic-localized transition for 
nonzero interactions. Additionally, our results suggest 
that the glassy-ergodic transition is a better candidate 
for a sharp phase transition. In spite of the fact that the 
local nature of the observable constitutes a limitation to 
perform a reliable finite size scaling, our strategy seems 
promising to analyze different underlying topologies and 
different ways to breakdown integrability. Last, but not 
least, given the state-of-the-art in Nuclear Magnetic Res- 
onance IS" 33 ; the evaluation of high temperature corre- 
lation functions, like the present version of the LE, could 
become experimentally accessible with the suitable elec- 
tion of the system. 
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FIG. 3. Phase diagram for A/i,i(r) at T = 12fi/J. Data point 
(a) is given for the MBL transition in Ref. 28], A = l.OJ, 
he — 2.7±0.3. Data point (b) is given for the MBL transition 
in Ref. [1], A = l.OJ, he = 3.5 ± 1.0. These points are 
slightly shifted in the plot from A — l.OJ in order to avoid 
their overlap. 
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